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On ' ^'^^Pb(^^Kr,n) reaction at Berkeley is systematically studied in the Relativis- 



The structure of superheavy elements newly discovered in the 



tic Mean Field (RMF) approach. It is shown that various usually employed 
RMF forces, which give fair description of normal stable nuclei, give quite 
different predictions for superheavy elements. Among the effective forces we 



^ ' tested, TMl is found to be the good candidate to describe superheavy ele- 

ments. The binding energies of the ^^^118 nucleus and its a— decay daughter 
nuclei obtained using TMl agree with those of FRDM within 2 MeV. Similar 
conclusion that TMl is the good interaction is also drawn from the calculated 
binding energies for Pb isotopes with the Relativistic Continuum Hartree Bo- 
goliubov (RCHB) theory. Using the pairing gaps obtained from RCHB, RMF 
calculations with pairing and deformation are carried out for the structure of 
superheavy elements. The binding energy, shape, single particle levels, and 
the Q values of the a— decay Qa are discussed, and it is shown that both 
pairing correlation and deformation are essential to properly understand the 
structure of superheavy elements. A good agreement is obtained with exper- 
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imental data on Q^^. 
PACS numbers : 21.60.Jz, 21.65.+f, 21.10.-k, 21.10.Gv, 27.90.+b 

Keywords: superheavy elements, Relativistic Mean Field, Pairing correlation, deforma- 
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I. INTRODUCTION 

Following the discovery of a— decay isotopes of Elements Z = 110, 111, and 112 at GSI 
an isotope of the Element 118, ^^'^118, and several its a— decay daughter nuclei were 
announced to have been discovered at Berkeley Lab's 88-Inch Cyclotron with the newly 
constructed Berkeley Gas-filled Separator by bombarding lead target with an intense beam 
of krypton ions of 449 MeV @] . The sequence of decay events is consistent with the long- 
standing theoretical prediction that there exists an "island of stability" around 114 protons 
and 184 neutrons and activates once again the study of superheavy elements. 

The study of superheavy elements has been a hot topic for the last two decades. Recent 
works on the collisions, structure and stability of Heavy and Superheavy Elements can be 
found in Refs. |5HTT|. In a recent paper, Smolanczuk claimed that the reaction "^^^Pbl^^Kr, n) 
should have a particularly favorable production rate |]12[. This motivated the experiment 
at Berkeley. According to the authors, the synthesized superheavy element ^^^118 decays 
by emitting an alpha particle within less than a millisecond, leaving behind the isotope of 
element 116 with mass number 289. This daughter nucleus is also radioactive, alpha-decaying 
to an isotope of element 114. The chain of successive alpha decays continues until element 
106. 

Smolanczuk discussed also the properties of superheavy elements in this mass region 
under the constraint of a spherical shape based on a macroscopic-microscopic approach [13 



In contrast to his approach, here we study the structure of superheavy element ^^'^118 and of 
the daughter nuclei in the sequence of a decays in the Relativistic Mean Field (RMF) theory. 
The effects of deformation and pairing correlation will be taken into account. The pairing 
gaps for deformed RMF calculations are taken from the Relativistic Continuum Hartree 
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Bogoliubov (RCHB) theory |T^, which is an extension of the Relativistic Mean Field and the 



Bogohubov transformation in the coordinate representation |jT5|. As the spin-orbit sphtting 



which governs the shell structure and magic number is naturally obtained in the RMF theory, 
we expect that the structure of superheavy elements can be understood properly once the 
deformation and pairing correlation are taken into account. We investigate the binding 
energy, deformation, the Q-values of the alpha decay, the effect of pairing correlation, shell 
structure, and the structure of single particle levels for protons and neutrons. 

The paper is organized as follows. In sect. II, we present the results of RMF calculations 
without pairing correlation for several standard forces, which give fair description of normal 
stable nuclei. We thus discuss the appropriate force to describe superheavy elements. In 
sect. Ill the RCHB theory is used to investigate the pairing correlation in these superheavy 
elements. The RCHB provides not only a unified description of mean field and pairing 
correlation but also a proper description for the continuum and the coupling between the 
bound state and the continuum |jT^. We then perform in sect. IV the study by a deformed 
RMF+BCS approach using the pairing gaps supplied by RCHB. We summarize the paper 
in sect.V. 

II. EXAMINATION OF VARIOUS RMF PARAMETER SETS 

There are many parameter sets for RMF calculations, which provide nearly equal quality 
of description for stable nuclei. Therefore, we wish to find at first which effective force in 
RMF is more suitable to describe superheavy elements. As claimed in Ref. the results are 



strongly interaction dependent. For this purpose we perform RMF calulations that include 
deformation but ignore pairing correlation with different effective forces. The details of the 



method can be found in Ref. 16 



Table | compares the binding energies E of superheavy element ^^^118 and its a decay 
daughter nuclei calculated with effective forces TMl NLl, NL3 and NLSH. For compar- 
ison, the results of the phenomenological FRDM calculations are given in the last column 
f|18i. 
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The results of TMl are nearly the same as those of FRDM for 265io4, 269^06, ^^'^108, 
^■^^110, and 281112. They are within 1 MeV from each other. The difference between TMl 
and FRDM results gets larger for 285ii4^ 289;^;^g ^^^^ ^^^n^^ jg g^j^ smaller than 3 MeV. 
Though there are differences of several MeV, NLl and NLSH give similar results as TMl. 
The NL3 parameter set, on the other hand, gives a difference of about 50 MeV from the 
other calculations. 

One important difference between the RMF calculations with TMl and FRDM is that 
the additional gain of the binding energy when one moves from 2^^112 to 285]^ jg much less 
in the RMF calculations. In other words, Z=114 has a weaker meaning as a magic number 
in the RMF calculations. Though, strictly speaking, it may not be adequate, let us call this 
effect as the change of the shell structure or of the magic number property at Z=112. As 
we see shortly, a similar effect appears in the Z dependence of the nuclear shape. This effect 
eventually plays an important role in reproducing the qualitative trend of the experimental 
data on the atomic number dependence of Qa- 

TABLE 1. The binding energy E of the superheavy element '^^^118 and of its a decay daughter 
nuclei calculated in the RMF theory with different effective interactions TMl, NLl, NL3, and 
NLSH. The prediction of FRDM is also given in the last column. 





TMl 


NLl 


NL3 


NLSH 


FRDM 


'^^104i6i 


1949.2 


1953.0 


1912.5 


1950.5 


1950.0 


269106163 


1970.2 


1974.1 


1932.0 


1977.5 


1970.5 


2^3108165 


1990.2 


1992.8 


1950.9 


1997.9 


1989.4 


2^^110167 


2007.9 


2010.5 


1968.2 


2016.1 


2007.0 


2«^112l69 


2026.2 


2029.8 


1986.2 


2034.7 


2025.2 


285ll4in 


2041.3 


2051.8 


2003.0 


2048.7 


2044.1 


2^^116173 


2058.2 


2068.7 


2014.2 


2065.6 


2061.1 


293118175 


2074.7 


2088.9 


2028.7 


2080.5 


2077.2 
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Similar trend concerning the mutual comparison of different forces appears also in Table 
0, where the Q-values of a decay sequence Q« = E{^He) + E{Z -2,N -2) - E{Z, N) 
(MeV) are shown. The Q-values given by TMl and FRDM are quite similar except for 
^^^114, where the difference is 3.8 MeV. This large difference is connected with the change 
of the shell structure mentioned above. 

TABLE II. The Q-values of a decay for superheavy element ^^^118 and its a decay daughter 
nuclei. We used the experimental value for the binding energy of the a particle. 



^Zr, 


TMl 


NLl 


iVL3 


NLSH 


FRDM 


^^^106163 


8.3 


7.2 


8.8 


1.3 


7.8 


^^^108165 


8.3 


9.6 


9.4 


7.9 


9.4 


2^^110167 


10.6 


10.6 


11.0 


10.1 


10.7 




10.0 


9.0 


10.3 


9.7 


10.1 




13.2 


6.3 


11.5 


14.3 


9.4 


^^^116173 


11.4 


11.4 


17.1 


11.4 


11.3 


^^^118175 


11.8 


8.1 


13.8 


13.4 


12.2 
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Table |T| shows the corresponding deformation parameter p in the ground state. TMl 
predicts a stable prolate deformation /3 ~ 0.2 for all the nuclei listed in the table, taking 
the minimum at ^'^^112. NL3 and NLSH give similar results as TMl, but the minimum 
deformation is shifted to ^^^114 for NL3. The NLl predicts a spherical shape for ^^^118, 
while FRDM almost spherical shape for ^'^'^110, ^^^112, ^^^lU, ^^Hl6 and "^^^IIS. The shift 
of the atomic number, where the deformation becomes minimum, from Z=114 to 112 is what 
we already mentioned as an evidence of the change of the shell structure. 

TABLE III. The deformation (3 of the superheavy element ^^^118 and its a decay daughter 
nuclei calculated with different effective interactions TMl, NLl, NL3, and NLSH. The prediction 
of FRDM is also given in the last column. 



Ay 


TMl 


NLl 


iVL3 


NLSH 


FRDM 


265l04i6l 


0.2656 


0.2687 


0.2685 


0.2513 


0.222 


^^^106x63 


0.2059 


0.2576 


0.2479 


0.2073 


0.221 


2^^108x65 


0.2021 


0.2438 


0.2032 


0.2068 


0.173 


2^^110167 


0.1857 


0.2279 


0.1865 


0.1914 


0.089 


'^^112169 


0.1681 


0.2072 


0.1683 


0.1741 


-.096 


285114m 


0.2118 


0.1593 


0.1606 


0.2316 


0.080 


289i 1 c 


0.2373 


0.1548 


0.2163 


0.2566 


0.081 


293i 1 Q 
110175 


0.2340 


0.0600 


0.2402 


0.2946 


0.080 
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Table |^ compares the corresponding charge-radii Re- In contrast to the big difference 
seen in the binding energy, the charge-radii Rc for different forces he within 1% from each 
other. 

TABLE IV. Comparison of the charge-radii Rc of superheavy element ^^^118 and its a decay 
daughter nuclei calculated with different parameter sets. 



Ay 


TA/f 1 

J. iVj 1 


NT ^ 


1 V ±jkj 


1\ ±JkJ 11 


26^106163 


6.146 


6.153 


6.177 


6.104 


2^^108165 


6.174 


6.176 


6.197 


6.134 


2^^110167 


6.197 


6.203 


6.220 


6.158 


2^^112169 


6.226 


6.228 


6.247 


6.188 


285ll4in 


6.282 


6.248 


6.275 


6.250 


2^^116173 


6.332 


6.272 


6.340 


6.303 


2^^118175 


6.360 


6.276 


6.392 


6.355 
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III. PAIRING CORRELATION IN SUPERHEAVY ELEMENTS: DESCRIPTION 

BY RCHB 



In this section we study the effects of pairing correlation in superheavy element ^^^118 
and its a decay daughter nuclei by using the self-consistent and fully microscopic RCHB 



theory [|T4[ under the constraint of a spherical shape. With the pairing gap obtained from 
RCHB, a self-consistent and more complete RMF calculation with both pairing correlation 
and deformation will be carried out in the next section. 

Before applying the RCHB theory to newly discovered superheavy elements, we examine 
once again which effective force is the most suitable to describe superheavy elements. For 
this purpose, we use lead isotopes as test cases. 

TABLE V. Comparison of the binding energies E of Pb isotopes calculated in RCHB theory 
with 4 different parameter sets with experimental data. The last four columns are the root mean 
square neutron, proton, matter and charge radii calculated with the TMl parameter set. 



A(Pb) 


Exp. [Me V] 


TMl 


NLl 


NL3 


NLSH 


Rn 


Rp 


Rm 


Rc 


202 


1592.20 


1592.91 


1596.60 


1565.43 


1596.00 


5.629 


5.420 


5.545 


5.479 


204 


1607.52 


1609.18 


1611.35 


1580.36 


1611.35 


5.656 


5.429 


5.566 


5.487 


206 


1623.40 


1623.78 


1625.73 


1594.77 


1626.11 


5.683 


5.437 


5.586 


5.495 


208 


1636.45 


1637.76 


1639.72 


1608.56 


1639.99 


5.713 


5.447 


5.609 


5.505 


210 


1645.57 


1646.92 


1646.78 


1618.04 


1648.34 


5.743 


5.467 


5.636 


5.525 


212 


1654.52 


1655.73 


1653.57 


1626.97 


1656.43 


5.772 


5.486 


5.663 


5.544 
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The binding energies of six Pb isotopes calculated by RCHB with 4 different effective 
forces are compared with experimental data in Table . Although all the calculations except 
for NL3 well reproduce the experimental binding energies of the Pb isotopes, TMl gives the 
best reproduction of the data. Therefore we expect that the RMF calculations with TMl 
and pairing correlation will give a satisfactory description of superheavy elements. The rms 
radii for neutron R^, proton Rp, matter Rm, and charge radii Rc calculated by RMF with 
TMl are given in the last four columns in Table |V|. 

TABLE VI. The binding energy E, one neutron separation energy Sn, the Q value for the a 
decay Qa, matter and charge rms radii Rm and Rc, neutron and proton pairing gaps in RCHB 
with TMl for the superheavy element ^^^118 and its a decay daughter nuclei. The results for ^^Kr, 
^''^Pb and ^^^118 are also given. 





E 




Qa 


Rm 


Rc 


A„/Ap 




750.2 






4.221 


4.182 


-0.010/-1.304 




1637.6 






5.649 


5.541 


-0.000/-0.000 


265l04i6l 


1944.4 






6.175 


6.099 


-0.622/-1.173 


2^^106163 


1965.6 


4.9 


7.1 


6.202 


6.131 


-0.421/-1.146 


'^^108165 


1986.8 


5.9 


7.1 


6.228 


6.160 


-0.283/-1.133 


2^^110167 


2007.2 


5.5 


7.9 


6.255 


6.189 


-0.380/-1.092 


2^^112169 


2027.1 


6.3 


8.4 


6.281 


6.218 


-0.338/-1.030 


285ll4i7l 


2046.8 


6.6 


8.6 


6.308 


6.247 


-0.030/-0.948 


289 1 1 c 
iibi73 


2065.2 


5.8 


9.9 


6.333 


6.273 


-0.013/-0.841 


293i 1 Q 


2082.1 


5.9 


11.5 


6.356 


6.296 


-0.315/-0.696 


ii»176 


2088.8 


6.7 




6.364 


6.299 


-0.442/-0.702 
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We have then calculated the binding energy E, one neutron separation energy Sn, the 
Q value for the a decay Qa, matter and charge rms radii Rm and Rc, neutron and proton 
pairing gaps for superheavy elements in the RCHB with TMl. The results are shown in 



Table The matter rms radius Rm is larger than the charge rms radius Rc for all nuclei 
due to the neutron excess. The Q value of the a decay increases monotonically with Z. The 
proton pairing gap parameter is around 1 MeV, while the neutron pairing gap parameter is 
relatively small due to blocking effects. The calculation for ^^Kr, ^osp^ and ^^^118 are also 
given for reference to understand the fusion barrier to synthesize the element ^^''llS. 

In Figs.|l| and 0, the single particle levels in the canonical basis for neutrons and protons 
in ^^^118 are given, respectively. In order to avoid the irregularity due to the blocking effect, 
we give the single particle levels in ^^^118 instead of ^^^118. The Fermi surface for neutrons 
and protons is given in each figure by the dashed line. The potential is the sum of the vector 
and scalar potentials. Fig.|l| indicates that, after the sub-closed shell at = 164, the next 
closed or sub-closed shells occur at = 198 and = 210. For the proton case, closed or 
sub-closed shells occur at Z = 106, Z = 114 and Z = 120. 

The Fermi level for protons in ^^'^118 is at A = —1.916 MeV, while that for neutrons at 
A = —6.304 MeV. Although the Fermi level for protons is very close to the continuum, the 
wave functions of all the protons are well localized in a small region because of the Coulomb 
barrier. 

Figs.^ and ^ show the change of the single particle neutron and proton levels near the 
Fermi surface along the a— decay chain from ^^^118. Similarly to Figs.|T] and 0, we give the 
single particle levels for the neighboring even-even nuclei in order to avoid the irregularity 
due to the blocking effect. Adding an a particle always raises the proton single particle 
levels and lowers the neutron single particle levels. There are distinct gaps of about 2 MeV 
at A^ = 164 and 172 and of about 3 MeV at A^ = 198 for neutrons. 

The a— decay energy is shown in Fig.^ as a function of the atomic number along the 
decay chain from ^^^118. The observed data and the prediction of FRDM are also included, 
where the former are taken from Fig.4 in |Q. Compared with the data, RCHB calculations 
give systematically too small Qa- This reflects the deformation effect which is ignored in 
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the present RCHB calculations. The Qa calculated in the RMF calculations neglecting the 
pairing correlation but including the deformation (the open circles) is somewhat larger than 
that in RCHB but still smaller than the data for Z < 108 and fluctuates for larger Z showing 
a sharp peak ai Z = 114. This contrasts to the result of FRDM, which also fluctuates, but 
shows a deep minimum at the same atomic number reflecting a sub-closed shell at Z = 114 
in this model. 




2 4 6 8 10 12 14 16 

R[fm] 

FIG. 1. The single particle levels in the canonical basis for neutrons in ^^^118 calculated by 
RCHB with TMl. The neutron potential Vy(r) + Vs{r) is represented by the solid line and the 
Fermi level by a dashed-line. 
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O <N CD 00 



[A9|A|]"3 

FIG. 3. The single partiele levels for neutrons near the Fermi surface calculated in the RCHB 
with TMl. They are shown for the neighboring even-even nuclei to the superheavy elements in the 
a— decay chain from ^^^118. 
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FIG. 5. Comparison of the theoretical a particle energies Qa with the observed data. 



IV. THE DESCRIPTION OF SUPERHEAVY ELEMENTS BY DRMF+BCS 

Using the pairing gap from RCHB, we now perform the RMF calculations by including 



both deformation and pairing correlation. The results are given in Table VII for the bindin; 
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energy E, the a— particle energy, matter and charge radii, and the neutron, proton and 
matter deformation parameters. The calculated binding energies for ^^Kr, "^^^Pb and ^^^118 
are also given. Each binding energy increases by 0.3 to 2 MeV with the pairing correlation, 
and can noticeably alter the atomic number dependence of Qa- 

We added the results of Qa calculated by the DRMF+BCS in Fig.|^. Comparing with 
RMF calculations, we observe that the theoretical Qa becomes much closer to the experi- 
mental data by the inclusion of the pairing correlation. Only for Z = 114, the Qa remains 
the same and has a difference of 2 MeV from the data. Interestingly, Qa takes maximum 
at Z = 114 in DRMF+BCS in accord with the hump in the data, while the FRDM gives a 
minimum there. 

TABLE VII. The binding energy E, a— decay Q value, matter and charge radii, and neutron, 
proton and matter deformation parameters calculated in the DRMF+BCS theory with TMl. 



""Zn 


E 


Qa 




Rc 


(3n 








751.0 




4.222 


4.189 


0.003 


0.005 


0.004 




1636.8 




5.650 


5.544 


0.000 


0.0000 


0.000 


265l04i6l 


1951.4 




6.206 


6.117 


0.212 


0.210 


0.211 


^^^106163 


1971.9 


7.8 


6.232 


6.147 


0.208 


0.207 


0.208 


'^^108165 


1991.1 


9.1 


6.256 


6.175 


0.197 


0.197 


0.197 


^^^110167 


2009.4 


10.0 


6.279 


6.202 


0.178 


0.178 


0.178 


^^^112169 


2027.3 


10.4 


6.302 


6.228 


0.164 


0.164 


0.164 


285ll4in 


2042.4 


13.2 


6.347 


6.280 


0.204 


0.202 


0.203 


iibi73 


2058.7 


12.0 


6.390 


6.328 


0.227 


0.228 


0.227 


293 n o 


2075.0 


12.0 


6.423 


6.363 


0.230 


0.234 


0.232 


294i 1 o 
ii»176 


2081.3 




6.429 


6.365 


0.224 


0.229 


0.226 
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V. SUMMARY 



We made a systematic study of the structure of superheavy elements recently discovered 
at Berkeley Lab's 88-Inch Cyclotron by the reaction ^^Kr + ^°^Pb at 449 MeV in the frame- 
work of Relativistic Mean Field (RMF) approach. We have shown that usually used various 
RMF forces, which provide fair description of normal stable nuclei, give quite different pre- 
dictions for superheavy elements. Among them TMl is found to be the good candidate to 
describe superheavy elements. 

We have shown that the binding energy obtained from TMl agrees with that of FRDM 
within a difference of 2 MeV. The same conclusion that TMl is the good interaction has 
been drawn from the calculations of the binding energy of Pb isotopes using Relativistic 
Continuum Hartree Bogoliubov (RCHB) theory. However, neither the deformation nor the 
pairing correlation alone could explain the data of Qa- 

We then performed RMF calculations of superheavy elements which include both the 
pairing correlation and deformation by using the pairing gaps obtained from RCHB. We 
have thus shown that a good agreement can be obtained between theory and experimental 
data concerning the Q value of the a-decay. Especially, our RMF calculations reproduce a 
peak at Z=114 seen in the experimental data. We conjecture that this peak appears because 
of the shift of the shell structure, e.g. concerning nuclear shape, from Z=114 in FRDM to 
Z=112. 

Finally we wish to make a few comments on open questions. We kept the pairing gap 
parameter once it has been fixed for a spherical shape, ignoring the possibility of its shape 
dependence [|T^. Another basic assumption is that the observed a— decays are from ground 
state to ground state, though this might not be the case for a part of the a— decay chain. 
We noticed a paper by Cwiok et al. after we have completed our study. The validity 
of the above mentioned approximation and assumption will be worth being examined to 
obtain more reliable understanding of superheavy elements. It would also be important to 
understand the difference between our conclusions and those in ref. [^, where the authors 
predict systematically much smaller deformation for all nuclei and also claim that the de- 
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formation monotonically decreases towards Z=118. We will address these questions in a 
separate paper. 
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